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In this letter, we show that the non-linearities of large amplitude motions in atomic nuclei induce 
giant quadrupole and monopole vibrations. As a consequence, the main source of anharmonicity is 
the coupling with configurations including one of these two giant resonances on top of any state. 
Two-phonon energies are often lowered by one or two MeV because of the large matrix elements 
with such three phonon configurations. These effects are studied in two nuclei, ''"Ca and ^'^*Pb. 
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Many-body fermionic systems possess collective vibra- 
tional states which are well described as bosonic modes 
(phonons). The existence in atomic nuclei of such states, 
both low-lying and Giant Resonances (Gil) is now well 
established up to the second quantum ||l|,||. However, 
their properties such as energy and excitation probabil- 
ity are still open questions. From the experimental point 
of view the strong excitation cross section of two phonon 
states calls for the presence of large anharmonicities but 
up to now, all the theoretical estimates were pointing 
to weak deviations from a harmonic spectrum. To our 
knowledge, so far only the mixing of one- and two-phonon 
states has been considered in microscopic calculations, 
with two exceptions. In ref. the coupling to some spe- 
cific three-phonon configurations has been included as a 
mechanism generating the damping width of the Double 
Giant Dipole Resonance. In ref. the fragmentation 
of the doubly excited low lying octupole states in ^°^Pb 
has been studied by allowing the coupling to one- and 
three-phonon configurations with a low energy cut-off in- 
troduced to reduce the diagonalization space. For this 
reason monopole, (GMR) and quadrupole (GQR) con- 
tributions which, as we will see, play an important role, 
were neglected. 

In the present paper we show that a correct description 
of the states for which the main component is a two- 
phonon configuration requires the inclusion of one- and 
three-phonon ones. We stress the essential role played by 
the breathing mode in the nuclear anharmonicity as an 
important novelty of the present analysis since volume 
modes are usually not considered in damping or coupling 
mechanisms. Moreover, we will show that the very col- 
lective GQR plays also an important role. 

The starting point of our calculation is a mapping of 
the fermion particle-hole operators a^a^ into boson op- 
erators as for example the one proposed in ref. 

Kh + (1 - E K'h'K'Ak' + .... (1) 

p' h' 

h p 



where (a) creates (annihilates) one nucleon in an oc- 
cupied (h) or unoccupied (p) single particle state. The 
second term on the right hand side of eq.(|l|) is a cor- 
rection that takes care of the Pauli principle. Then we 
construct a boson image of the Hamiltonian, truncated at 
the fourth order in the and B operators. Introducing 
the Bogoliubov transformation for bosons: 

Qt ^ J2(^phKh - VhBph) (3) 

and imposing that the quadratic part of the boson Hamil- 
tonian in the new operators is diagonal, we obtain the 
usual Random Phase Approximation (RPA) equations 
for the X and Y amplitudes. 

By inverting eq.(|^) we can express Hb in terms of the 
collective and Q operators: 

Hb = HnQ^Q + (HziQ^Q^Q + h.c.) + H22Q^Q^QQ 
+ (H3og^g^Q^ + h.c.) + {nsiQ^Q^Q^Q + h.c.) 
+{^AoQ^Q^Q^Q^ + h.c.) + ... (4) 

with Till = E^6^,y'. The Ti matrices are expressed in 
terms of the X and Y of transformation (^). The con- 
tributions to eq.(^) coming from the high order terms of 
the expansion (^) appear to be reduced by the number of 
configurations involved in the collective states and there- 
fore can be neglected. In the case of closed shell nuclei, 
the RPA correlations are moderate. Therefore, the Y/X 
ratios are small. In eq.(||) we will neglect the Ti terms 
containing at least one Y amplitude. These two approx- 
imations leave unaffected only the first three terms of 
eq.(^) 1^. We will compare the spectra of ^"Ca and ^"^Pb 
obtained by the diagonalization in the spaces containing 
up to two-phonon states and up to three-phonon states, 
respectively. In ref. |^ a similar analysis was done in the 
two level Lipkin model. It was found that this approxi- 
mation is well justified and one gets good results in the 
larger space for the eigenstates which main component 
is a two-phonon configuration. 

All calculations have been performed by using the SGII 
Skyrme interaction We include all natural parity 
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RPA collective one-phonon states with angular momen- 
tum J < 3 which exhaust at least 5% of the EWSR 
and all two- and three-phonon configurations built with 
them, without any energy cut-off, with both natural and 
unnatural parity. 

Let us start looking at the results for ''"Ca. In table || 
we show the one-phonon states taken into account. In 
table |l| we show some results of the diagonalization for a 
selected set of states. The energies obtained in the space 
up to two-phonons (see ref. 1^) are reported here for 
comparison. The so-calculated anharmonicity was lim- 
ited to a few hundred keV. Let us now study the more 
complete calculation including the three phonon states. 

As a general comment, one can say that the shift in- 
duced by the coupling to three-phonon states is fairly 
large, being in almost all the cases more than 1 MeV, 
and always downward. This can be understood in second 
order perturbation which, as can be seen from the table, 
gives a good estimate for the energies in most cases. In 
second order perturbation the correction to the energy is 
given by 



AS, ipi\V\ipi > 



EO 



(5) 



where \ipi > is the considered unperturbed state, \ipj > 
all the other states and the corresponding unper- 
turbed energies. Since the diagonal, first order, contri- 
bution is small in most cases, the sign of the shift is that 
of the denominator in the second order term. Therefore, 
if > is a two-phonon state, the contributions from 
three-phonon configurations are negative in most cases 
since most of the three phonon states lye above the two 
phonon ones. Moreover, whenever a GMR is added on 
top of any state, the H21 terms (see eq.^) are large , of 
the order of 1 to 2 MeV in *'^Ca. The specific values 
can be found in the last three columns of table |. Indeed, 
in the coupling leading to the addition of one GMR on 
top of any state, the residual interaction between the un- 
derlying fermions is not truncated by conservation laws, 
because the particles and the holes involved in the GMR 
carry identical parity and spin quantum numbers (cf. ref. 
1^). Phenomenologically, this strong coupling of all col- 
lective vibrations with the breathing mode comes from 
the fact that in a small nucleus like the ''°Ca any large 
amplitude motion affects the central density. Therefore, 
surface modes cannot be decoupled from a density varia- 
tion in the whole volume as clearly seen in recent TDHF 
simulations ref. pT] . 

If the state is a two-phonon one, then the matrix el- 
ements coupling it to the state obtained by exciting a 
breathing mode on top of it are about 3 MeV (up to 5.5 
MeV) when the less (more) collective component of the 
^"Ca GMR is considered. Even larger matrix elements 
are obtained, when the states connected by 7^21 involve 
several GMR. In that case a Bose enhancement factor 



appears and no Clebsch-Gordan coefficients enter in the 
calculation. Thus the matrix element between the dou- 
ble and the triple GMR located at 18.25 MeV, Mi, is 
\/6 times larger than between the single and the double 
Ml. That gives a matrix element of -5.22 MeV, giving 
a contribution of -1.49 MeV to the second order energy 
correction of the double Mi state. An even larger value 
comes out in the case of the double GMR located at 
22.47 MeV, A/2, and the triple M2, namely a matrix el- 
ement of -9.69 MeV giving a -4.18 MeV contribution to 
the energy shift of the double M2. This is due to the fact 
that M2 is more collective than Mi in ^"Ca. 

Something similar, but less strong, happens also for 
the matrix elements connecting some state with that 
built by adding one GQR phonon. We quote two ex- 
amples. The low-lying component of the Giant Dipole 
Resonance \Di > has a matrix element of the residual in- 
teraction with the states \Di Mi >, \Di (g) M2 > and 
\Di (g> Qi > equal to -1.38 MeV, -2.12 MeV and -1.25 
MeV respectively. Another example, with total J— I, 
is given by the matrix elements between jZ?! Qi > 
and I (Ml Di)i <S) Qi >, |(M2 «) Di)i (E) Qi > and 
\{Qi)l ® Di> equal to -2.74 MeV, -4.61 MeV and -1.41 
MeV respectively. 

These findings clearly indicate that large amplitude 
motions are strongly coupled both to surface and volume 
oscillations, the latter being more important in ^"Ca. It is 
worthwhile stressing that such large corrections to the en- 
ergy of two-phonon states are obtained despite the quite 
large absolute values of the energy difference between the 
coupled states. Therefore, introducing an energy cut-off 
in the three-phonon states included in the calculation 
may lead to erroneous results. Let us consider for exam- 
ple the case of the 0+ member of the multiplet of double 
low-lying octupole states. At first order perturbation, 
it is shifted up by 2.24 MeV. The second order correc- 
tion coming from the single GMR states is -0.93 MeV. 
These two contributions, leading to a total shift of +1.31 
MeV, dominate the effects of the coupling with one- and 
two-phonon states as confirmed by the result of the diag- 
onalization in this subspace. When three-phonon states 
are included, one gets a further shift down of 1.86 MeV 
coming from the configuration including a GMR on top 
of the two octupoles. This contribution is absent in ref. 
Q because the energy cut-off introduced there in order 
to reduce the number of three-phonon configurations was 
too low. The same happens for the other members of the 
multiplet as well as for the double Di or I?2, the double 
Qi and the I?ior D2 ® Qi states. 

The results for 208pj^ shown in tables III and |rv| . 
The same general remarks already made for '^'^Ca apply 
also in this case. The most relevant difference is that the 
role played by the GMR and the GQR in ^"Ca IS now m- 
verted, the latter being dominant in ■^°*Pb. This reduced 
importance of the GMR may come from the fact that 
in large nuclei the surface vibrations can occur without 
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changing the volume. Concluding about the energy of the 
two-phonon states one can see that the inclusion of the 
three phonon configurations induces an anharmonicity of 
more than 1 MeV in *°Ca but only of a few hundred keV 
in ^°*Pb. Because of the location at high energy of the 
three phonon states, the observed shift is systematically 
downward. It is important to stress that the considered 
residual interaction only couples states with a number of 
phonon varying at maximum by one unit. Therefore, the 
energy variation of the two-phonon spectrum induced by 
inclusion of four and more phonon states would be small 
since it corresponds to a third order perturbation involv- 
ing two large energy differences in the denominator. 

If we now analyze the splitting of the two-phonon mul- 
tiplets we can see that it remains small for giant reso- 
nances (about a few hundred keV) while it may go up to 
1 MeV for low lying states in ^"Ca. Comparing the split- 
ting and the ordering of the states obtained in first order 
perturbation and in the full calculation we can see that 
they remain almost unchanged. Therefore, the diagonal 
matrix elements of the residual interaction are responsi- 
ble for this splitting and ordering. 

Let us now investigate the mixing induced by the resid- 
ual interaction. In tables |l| and [V the mixing coefficients 
of the two main components in each state are presented. 
First we can see that there is always one component that 
remains very large, explaining the success of the pertur- 
bation approach. The important point is that in general 
we observe large mixing coefficients, namely about 0.2 
to 0.4 or more in ^OCa and 0.15 to 0.3 in ^ospb. This 
may have very important consequences in the excitation 
process as we will investigate in a forthcoming work. 

It is worthwhile mentioning that, in some cases, a 
three-phonon component appears with a large amplitude 
in the wave function of a (mainly) two-phonon state, de- 
spite the fact that the residual interaction does not couple 
directly these configurations together. In a few cases, in- 
deed, this is the second main component as can be seen 
for 4'^Ca in table || (the \ (Di)l > and |M2 Qi > states) 
and for ^ospb in table |l^ (the \{Mi)^ > state). This 



happens because the diagonal matrix elements of the 
Hamiltonian in the two-phonon and three-phonon con- 
figurations are close and the matrix elements coupling 
the latter with other configurations are large. A similar 
situation has been found in ^"^^Pb for two one- phonon 
(mainly) states which have a three-phonon configuration 
as second important component, even though our Hamil- 
tonian does not couple directly states which numbers of 
phonons differ by more than one. This is the case for the 
state which main component is \Mi >, with amplitude 
Co= —0.79, and for which the second most important 
component is |(3^)2 ® 2+ > with ci— 0.55. Both com- 
ponents have large matrix elements with the two-phonon 
state |(3^)o >. How this mixing of the monopole res- 
onance may affect the monopole response, and so the 
usual conclusion about the compressibility, is now un- 



der study. The other case is the single high energy oc- 
tupole resonance \0 > which is strongly mixed with the 
states 1(2+ (8)3^)j i^Qi >. The energy of these states 
are, however, shifted by less than 100 keV. This is co- 
herent because the strong mixing is coming from a quasi 
degeneracy of the considered states. 

Summarizing, the spectrum of two-phonon states is 
strongly modified by their coupling to the three-phonon 
ones. All of the states appear mixed with the excitation 
of a GMR and GQR on top of it. This is due to the 
fact that most of the matrix elements of 7^21 coupling a 
phonon with the same phonon plus a GMR or a GQR are 
large. Moreover, because of the Bose enhancement fac- 
tors, the effect of 7^21 between the two and three phonon 
states is even larger. It is also to be noted that many 
of the important three-phonon states are higher in en- 
ergy than the two-phonon ones. Therefore, they induce 
a systematic shift down of the two phonon states as the 
sum of several quite large negative contributions. This 
unexpected finding can be understood as a modification 
of the central density in large amplitude motion leading 
to an excitation of the breathing mode. The case of the 
GQR seems to be related to the extreme collectivity of 
this state leading to a strong quadrupole response to the 
quadrupole component of the non-linearities of the mean- 
field. We also want to stress that, because of the pertur- 
bative nature of the observed phenomenon, the possible 
introduction of four-phonon states should not modify the 
above conclusions about two phonon states. Of course, 
our findings imply that in order to get a correct three- 
phonon spectrum one should further enlarge the space 
up to four-phonons. This is a formidable task which is 
beyond the scopes of the present paper. 
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TABLE III. Same as table I for the nucleus ^°*Pb. 



TABLE I. RPA one-phonon basis for Ca. For each state, 

spin, parity, isospin, energy and percentage of the EWSR are 
reported. In the following columns, Vmi stands for the matrix 
element < vlVlf ® Mi >, where v is the one phonon in the 
1"* column, the same for Vmj and Vq^. 
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O2 
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31.33 
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-0.31 



TABLE II. Results for Ca. In the first column, the states 
are labelled by their main component in the eigenvector and 
their unperturbed energy (in parentheses). In the second col- 
umn, the amplitude of the main component cq. Then for each 
total angular momentum J, we show the results of the calcu- 
lation in the basis up to 2 phonon states, the present results 
for the basis extended to 3 phonon states, the corresponding 
first order perturbation theory energy, and the second order 
one. The last two columns contain the 2"'' main component 
in the eigenstates and the corresponding amplitude ci. The 
subindex in the two-phonon configurations denotes J. All en- 
ergies are given in MeV. 
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TABLE IV. Same as table II for the ^°*Pb nucleus. 
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